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Abstract 

Some bright cD galaxies in cluster cooling flows have Thomson optical depths exceeding 
10~ 2 . A few percent of their luminosity is scattered and appears as diffuse polarized emission. 
We calculate the scattering process for different geometric combinations of luminosity sources 
and scattering media. We apply our results to galaxies, with and without active nuclei, 
immersed in cooling flows. We model observations of NGC 1275 and M87 (without active 
nuclei) in the presence of sky and galactic background fluxes which hinder the measurement 
of the scattered light at optical wavelengths. Current instruments are unable to detect the 
scattered light from such objects. However, when a galaxy has an active nucleus of roughly 
the same luminosity as the remainder of the galaxy in V, both the total and polarized 
scattered intensity should observable on large scales (5-30 kpc), meaning intensity levels 
greater than 1% of the background level. For typical AGN and galaxy spectral distributions, 
the scattering is most easily detected at short (U) wavelengths. We point out that a number 
of such cases will occur. We show that the radiation pattern from the central nuclear region 
can be mapped using the scattering. We also show that the scattered light can be used to 
measure inhomogeneities in the cooling flow. 
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I. Introduction 

A variety of observations suggest that large quantities of plasma and dark matter form 
halos around galaxies. Continuum X-ray brightness profiles (Jones and Forman 1984) demon- 
strate the existence of hot, distributed gas while X-ray line emission (Stewart, Canizares, 
Fabian, and Nulsen 1984) and centrally peaked X-ray images (Canizares, Stewart, and Fabian 
1983) show that cooling of the gas occurs near the center of some galaxies (see Fabian, Nulsen, 
and Canizares 1991 for a review of cooling flows). Optical observations of rich clusters show 
substantial line emission emanating from filaments of gas in the cluster core and/or from 
the central dominant galaxy (Kent and Sargent 1979, Heckman 1981, Cowie, Hu, Jenkins, 
and York 1983). Deprojection and modeling of the X-ray data (Stewart, Fabian, Jones, and 
Forman 1984 and others; a review is presented by Sarazin 1986) yield density profiles, mass 
inflow rates, and temperature distributions. Some of the mass accretion rates are enormous, 
e.g. 1000 M /yr for the cooling flow of PKS 0745-191 (Fabian et al. 1985). 

The inferred electron distribution implies Thomson optical depths of 10 3 10~ 2 over 
characteristic length scales (impact parameters) ranging from a galaxy radius of ~ 20 kpc 
(or inner edge of a cooling flow) up to a cluster core radius of ~ 100 kpc (or outer edge of 
a cooling flow). The galaxy's scattered luminosity is roughly 0.1-1% of its total luminosity 
and may be observable in some circumstances (Syunyaev 1982; Gil'fanov, Syunyaev, and 
Churazov 1987). Some galaxies with cooling flows have active nuclei, e.g. for a selected 
sample of eight radio-loud quasars, Crawford and Fabian (1989) find that all eight reside in 
clusters surrounded by hot ionized gas. Study of the underlying light distributions from radio- 
loud quasars also indicates that many are giant elliptical galaxies (Hutchings, Crampton, and 
Campbell 1984) for which cooling flows are expected. In general, massive galaxies with active 
nuclei are prime prospects for observing scattered light as suggested by Fabian (1989). He 
pointed out that a dense cooling flow around a quasar might yield an anisotropic optical 
continuum by the following mechanism. If the quasar's optical emission is beamed and 
collinear with its radio emission, then it scatters in the hot plasma of the cooling flow and 
produces continuum emission aligned with the radio axis, as has been observed (Djorgovski et 
al. 1987; Chambers, Miley, and van Breugel 1987; McNamara and O'Connell 1993). Strong, 
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frequency independent polarization of the scattered radiation is the characteristic signature 
of the process. In a study of six high redshift galaxies di Serego Alighier, Cimatti, and 
Fosbury (1992) found four show a high degree of linear polarization. Because the alignment 
of this polarization was perpendicular to the radio/optical axis they concluded that the 
polarization is caused by light, from an anisotropic active nucleus, scattered into the line of 
sight by gas surrounding the galaxy. 

For general wavelengths, Wise and Sarazin (1990) conclude that scattering is most likely 
detectable in the radio. If the scattering is measurable, Wise and Sarazin show that it can be 
used as a cosmological distance indicator as well as a probe of the physical conditions of the 
cooling flow. Recently, Wise and Sarazin (1992) examined the effects of inhomogeneities in 
cooling flows and temporal variations in the luminosity of the central source. The scattered 
radiation can provide a record of the luminosity over the last 10 5 years and may help identify 
the parent population of high-luminosity active galactic nuclei. 

The drawback with observing scattering at radio wavelengths is that Faraday rotation 
effectively erases the polarization. Distinguishing the scattered light from other diffuse emis- 
sion processes may be problematic once the polarization, the characteristic signature of the 
scattering process, is lost. At optical wavelengths Faraday rotation should be much less of 
a problem but detection may be more difficult. In this paper we investigate the intensity 
and polarization properties of visible light scattered by the typical, inferred electron column 
densities. In our model a bright galaxy (with or without an active nucleus) is surrounded 
by a halo of ionized matter, either a cooling flow or hot gas confined by a cluster potential. 
We calculate the intensity of the scattered emission, the polarization fraction, and the two- 
dimensional pattern of polarization seen on the sky. While some polarized light is always 
produced in these situations, we discuss the observational difficulties associated with its de- 
tection. We conclude that scattered, polarized light at optical frequencies should be detectable 
on large scales (5-30 kpc) in galaxies with sufficiently luminous nuclei and sufficiently dense 
cooling flows. The rough criteria are that the apparent luminosity of the nucleus exceed that 
of the parent galaxy and that the Thomson optical depth be of order 1 %. 

Observation of scattered light probes (i) the angular distribution of the radiation from 
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the nucleus and (ii) the angular and radial distribution of material that surrounds the il- 
luminating source. The first issue is especially relevant to theories that attempt to unify 
disparate classes of active galaxies with geometric beaming arguments (Barthel 1989). The 
second issue is important because the scattered polarized light can be used to probe length 
scales in cooling flows that have hitherto been inaccessible. X-ray observations have been 
the traditional tool for studying cooling flows, and much of our knowledge of them comes 
from the spectral and angular distribution recorded by the Einstein Observatory and more 
recently by ROSAT (Sarazin, O'Connell, and McNamara 1992). Here we show that scat- 
tered light provides a complementary source of information at angular scales smaller than 
the resolution of X-ray telescopes. 

The plan of the paper is as follows. In §IIa-f we review the physics of scattering. We 
lay out our assumptions and approximations and solve several simple, model problems. In 
§IIIa we show that the finite size of the luminosity source limits the fractional polarization 
at small distances from the center of the source. In §IIIb-f we carry out calculations for the 
prototypical sources M87 and NGC 1275, cD's with active nuclei, and an atypical source, 
PKS 0745-191, and discuss limitations imposed by diffuse emission from the parent galaxy 
and the sky background. Readers primarily interested in applications will want to scan §IIa 
and §IIf and then continue with §111. In §IV we summarize our work and future possibilities 
for research. 

II. Physics of Scattering 

a. Formalism 

We begin by assuming that a point source emits completely unpolarized light. Later, 
we will consider sources of finite extent. The scattering is assumed to be solely due to non- 
relativistic electrons. Although heavy elements will be present if the gas is not primordial, 
the characteristic temperatures of cooling flows are large enough that grains will not sur- 
vive except in the regions where star formation has begun. We ignore a photon's typical 
frequency shift bvjv ~ (A;T/m e c 2 ) 1//2 = 0.13Tg^ 2 , so that the radiative transfer is frequency 
independent. We also assume that the optical depth is everywhere small, so that single- 
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photon scattering is a good approximation. Finally, we assume that the luminosity of the 
central source is constant in time. 

The basic geometry is illustrated in Figure 1. Our discussion follows the standard treat- 
ment of Thomson scattering (e.g. Rybicki and Lightman 1979). Let the source be at the 
origin. Let the direction of the light emitted be k and let the direction of the scattered light 
be n. The plane in which k and n lie is the plane determined by the position of the source, 
the scattering material, and the observer. We define the first polarization vector, e\, as the 
direction perpendicular to k in that plane (or kx (nxk)) and the second, £2, as the direction 
perpendicular to h and k (or k x h). 

The differential Thomson cross section is 

Ir - M ~ ^ ' "> 2 > ■ « 

where r e is the classical radius of the electron. The unpolarized light emitted by the source 
may be decomposed into two components of equal intensity along e\ and €2- The scattered 
light is polarized because e\ • n is, in general, not equal to £2 • n = 0. Let 6 be the angle 
between k and n. The observed light, that is the light scattered into the line of sight, is de- 
scribed by two polarization directions, each perpendicular to n. We refer to the components 
perpendicular and parallel to k by _L and ||, respectively. The "emissivities" for each mode 
at position s are 

(3 = (I) fine 00s 2 9, (2) 

where F is the flux from the source and n e is the electron density. 

The intensity seen by the observer is the integral over the line of sight (I) of the scattered 
emission in each polarization mode. The two polarization directions are defined by 

ej_ = e\ x n 

e 11 = e± x n . (3) 

(In the quadrupole tensor, the signs of ej_ and are irrelevant.) Figure 2 illustrates the 
polarization directions as viewed by the observer. They lie along the "tangential" (ej_) and 
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"radial" (ey) directions of a circle whose center coincides with the luminosity source. They 
are the same at each point along a fixed line of sight. The intensities are 

3La T r°° n e (s) 



t r ■ 6LaT r 



oo s 2 



-dl 



/2 



r°° . 3La T /■«> n e (s)l 

,-J^ )dl = 64^7-00 — 

where L is the luminosity of the source. The fraction of polarized light is 

IjL-h 

n 



(5) 



|/± + /|| 

For a single source the observed plane of polarization lies along e±. 

h. Point Source with Centered Power-Law Gas Distribution 
Now we consider a very simple scattering situation, a power-law electron distribution 

n e (r) = n c (^j , (6) 

centered on the luminosity source. The total intensity of the scattered light at impact 
parameter h is 

T _ T n c a T /r e \" 3(a + 3)r(^) 

h \h) i287r 3 / 2 r(^) ' 1 ' 

Here, T is the Gamma function. The intensity of the scattered light is of course dependent on 
the flux of the source and the optical depth at the distance h. This expression for the total 
scattered intensity is identical to eq. (2.2) of Wise and Sarazin (1992) modulo a difference 
in units of 4ir. We see that / ~ F(h)r(h) where F(h) ~ flux and r(h) is the optical depth 
~ (n c crTh)(r c /h) a . The fraction of polarized light is 

n = — = — ^ — = ^—^ . (8) 

r n 2 /2 (l + sm 2 (9)) cos 2 +*(9)d9 3 + a 

For a broad range of power laws a significant fraction of the scattered light is polarized. 

The observer measures the Stokes parameters relative to some sky coordinate system (x, 
y) (see Figure 2) of his choice. From our results above, the scattering from a point source is 

/ = /x + /,l 
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Q = (7j.-/||)cos2* 

U = (I ± -7||)sin2* 

V = , (9) 

where \I/ is the angle between £ and ej_. For multiple sources of quasi- monochromatic 
radiation (I_u, I\\., and the Stokes parameters add up to give (Itot, Qtot, Utot, Vtot)- The 
total intensity seen at a point is Itot, the total polarization is 



Hot 

and the angle of polarization is given by 



T-r ^ QtOt + Uf 0t + Vf Qt 

v hot = j , l-LUJ 



tan2tt = ^. (11) 

Qtot 



c. Point Source with OS-Center Gas Distribution 



i.) Analytic Limiting Forms 

In general, the gas responsible for the scattering is hardly likely to be symmetrically dis- 
tributed about the illuminating source. In this section we relax the assumption of symmetry. 
We illustrate how displacing the luminosity center from the gas center alters the polarization 
fraction, the intensity of the scattered light, and the geometry of the polarization pattern. 
Assume an electron distribution with a scale-free power-law profile (c/. eq. (6)) with a = 1) 
centered at the point fb = (xo, yo, zq) and a luminosity point source at the origin. Without 
loss of generality, we assume that the observer views the scattered light along the z direction. 
We derive limiting forms for the scattering and polarization for lines of sight passing near 
fo and the origin, where the scattered intensity has local maxima. Elsewhere, the paucity of 
electrons and/or the dilution of flux from the source decreases the scattered intensity. We 
provide a systematic expansion in angular harmonics that describes the asymptotic behavior 
far from the center. 

Let the impact parameter of the line of sight with respect to the origin be p = \J x 2 + y 2 
and with respect to the density peak be q = \J (x — xq) 2 + (y — yo) 2 . Letting Iq = 
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3Ln c a T r c /6An 2 , eq. (4) becomes 

© = h II dl (X)wTPwhr^m- (12) 

First, for lines of sight that pass near the density peak (q/zo — > and finite p/zo), we expand 
in powers of (q/zo) and integrate term by term to find 

lL _ ( _ 21n ( ±) + ln4 u + r ^ _ 2tai r'fe) + o 





(13) 

Formally, ij_ and J|| diverge logarithmically near the density peak, although one (or both) 
of the assumptions, optically thin radiative transfer and unbounded growth in density near 
fo, must break down. As q — > the fractional polarization approaches the limit 







When the peak and luminosity source lie in the same plane (p/zq — > oo), the fractional polar- 
ization is a maximum. At a fixed-impact parameter, a larger line-of-sight separation lowers 
n. To summarize, significant fractional polarization is expected near a density peak whose 
three-dimensional displacement from the luminosity source is less than a modest multiple of 
the impact parameter. 

A systematic expansion for scattering at large impact parameters may be developed as 
follows. We express the density field (centered at tq with the form l/\r — ro\) in spherical 
harmonics and powers of the distance from the origin. The integral over the line of sight 
is converted to an angular integral and carried out term by term. Let (ro, 9o, 0o) be the 
spherical coordinates of the density field center, measured with respect to the origin (the 
luminosity point source). The impact parameter of the line of sight in the z, direction is p at 
angle 0. For p > rg we have the exact solution 



M A. E ( / H''"' : " "' : 7y(<WA,)-^. (15) 

^ImJ P 
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where 

Aim-1{ 1J (2 / + i)!! 



Jo £> 2 3p 3 TKJ ' ' 5p 4 



Sz m - ^1 ^7^) J 4 m . (16) 

To order I = 2 these expansions yield the following explicit forms: 

sin fl cos(0 - O ) + ^ [sin 2 6 cos(2(0 - O )) - - (3 cos 2 6 - l) J 

To = 3^ + l^ Sin ^° C ° S(0 ~ 0o) + 3^ (sin 2 ^ocos(2(0-0o)) + \ (3cos 2 £ - l)) ■ (17) 
The fractional polarization at large distances (to order / = 1) is 

- + ^sinfl o cos(0-0 o ) • (18) 

To summarize, at large distances from the luminosity source and density peak, the scattering 
approaches the simple form for a point source with a centered power law (cf. §IIb). As 
p — > 00, n agrees with eq. (8) for the adopted density power law (a = 1). 

The term by term integration may also be carried out for lines of sight that pass within 
the spherical shell of radius p < vq, but in this case there are three separate regions of 
integration and the formulae are more complicated. A consistent expansion is most easily 
derived using cylindrical coordinates. To summarize the steps: expand 1 / 1 r — r 1 using 
modified Bessel functions (I m and K m ), perform the line-of-sight integral in z, and then 
substitute a series expansion for I m (8.445 of Gradshteyn and Ryzhik 1980; hereafter G&R) 
and an integral definition for K m (G&R 8.432.3). An asymptotic expansion of the integrand 
gives integrals which may be reduced term by term to Beta functions. The results can be 
compactly expressed using s = ^ (x 2 + y 2 ) / (xq + yf) and t = ^Jzq/(xq + yfi); we assume the 
ordering s, t < 1. This ordering means that (i) the line of sight passes much closer to the 
luminosity source than the density peak and (ii) the density peak's displacement along the 
line of sight is less than or equal to its projected separation (other orderings may also be 
analogously treated). Let 

' s m+2k+l+8 \ 00 (_2) n 



T(m + k + U + 1 2 + 6 )T(k + H+1 2 +6 ) ((s - it) n + (s + it) n ) (19) 
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and let dj = 1 for i > and oq = 1/2: 

Ui = ?S> O " " m( *"* ,) l fi ^ tt -fl.(«.M)J- <20) 

Using the above expressions it is straightforward to give Iu and ij_ to order s /is t^ t for 
Us + Ht < M for a given M. For M = 2 the result is independent of t, 

1 1 S7T / 2 \ 

-T = -o l--s + scos(0-0 o ) (21) 

and 

i = ^(l-i S + S cos(0-0 o )). (22) 

The divergence in intensity as p — > may be traced to the assumption of a luminous point 
source. Note that the limiting behavior of the fractional polarization is 

n^i(l + ^-2 S cos(0-0 o )), (23) 

which is the same as for a point source embedded in a scattering gas of uniform density. To 
summarize, the fractional polarization expected in the vicinity of a luminosity point source 
approaches that of a point source in a uniform density medium. 

ii.) Numerical Results 

We now proceed with a more detailed description via a numerical calculation. We have 
carried out the quadratures in eqs. (4) for a grid of 100 by 100 lines of sight. As in the analytic 
work above we take a scale-free power law (eq. (6) with a = 1) for the density distribution. 
The source is located in the plane of the gas center (co-planar geometry zq = 0) or displaced 
out of the plane (\Jxq + y^/ zq = 1/4). Figure 3a shows the polarization fraction (length 
of the line segments), the angle of polarization (direction of the line segments), and the 
intensity of the scattered light (contours) for the co-planar geometry. Figure 3a and future 
figures are centered about the gas-density maximum (fo). The maximum intensity occurs 
near the source while a local maximum occurs near the gas center. The numerically derived 
fractional polarization, shown in Figure 3b, is a minimum near the source (~0.35; compare 
eq. (23)) and a maximum near the gas center (~0.75; compare eq. (14)). In the vicinity of 
the gas center and the point source, the polarization ranges between these extremes; at large 
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radii it approaches 0.5, as would be expected for a point source embedded in a gas following 
r _1 power law (eqs. (8) or (18)). Of course, at the gas center the numerically determined 
fractional polarization should approach its maximum value of unity, while at the source it 
should go to 1/3. These limits are only approximately observed in our numerical model 
because of the finite resolution in the 100 by 100 grid. When the source and gas center are 
displaced out of the plane (Figure 3c) scattering near the gas center is diminished (the flux 
from the source at the gas center is less; compare eq. (13)) and the polarization decreases 
(the scattering angle, 9, becomes small at the high-density region; compare eq. (14)). 

d. Comparison of Center and OS-Center Gas Distributions 

We now consider the scattering when both a centered and off-centered distribution are 
present. We want to derive simple criterion to distinguish which type of scattering domi- 
nates. Of course, since the radiative transfer is linear in our approximation the individual 
results may be simply added to give Itot and Utot as described in §IIb. First, assume both 
distributions lie in the same plane and that the off-center distribution is a power law with 
a — 1. The intensity from the off-center "clump" is given by 



jd _ J i.. 

1 (l- M 2 )V2 p 2 



1 + Vl-u 2 



1 - Vl -u' 2 



;!_ M 2)3/2 p 2 y 2 



1 + Vl-u 2 



(24) 



' 2 . 

where u = q/p, and q, p, and Iq have the same definitions as above. In the limit that the 
clump is centered on the galaxy, u — > 1 results in eqs. (7) and (8) are retrieved. 

For comparison, the ratio of I d = lf_ + I d , the total scattered intensity from the off- 
center clump, to I 9al , the intensity from an assumed galactic halo given by eq. (7), at the 
same impact parameter, is 

I d _ 3f(u) ( nfr d \ 



jgal g y n gal r gal J 

where f(u) is a function which diverges logarithmically at small u (this ignores the core 
radius of the clump). n g c al , n d , r gal , and r d are the scaling densities and radii used in eq. 
(6) for the off-center clump and the galaxy's gas distribution at the center. For fixed column 
densities of the clump and galaxy, the ratio depends only upon f(u). At a fixed distance from 
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the clump center (e.g. its core radius), u — > as the clump displacement jfrom the galaxy 
center increases. The result suggests that the clump scattering dominates halo scattering at 
large galactic radii but is sensitive to the large extent assumed for the clump (a — 1). 

It is straightforward to evaluate I cl / I gal for arbitrary power laws for galaxy (a gal > 0) 
and clump (a d > 2) when the line-of-sight contribution to scattering is dominated by the 
points of closest approach. We find the intensity ratio scales as 

I cl ^ ( n^rf 

Tqal ~ \ qal qal 

For the massive cD galaxies with which we are concerned, a gal ~ 1. Thus, for a sufficiently 
concentrated clump profile (a cl > 2), we see that at fixed distance from the clump center, the 
ratio of column densities fixes the ratio of the brightnesses. To summarize, our comparison 
shows that the most significant factor dictating the clump-to-halo brightness ratio is just 
the column density ratio. If the off-center clump has an extended (a = 1) distribution, 
the brightness ratio increases slowly (logarithmically) as it is displaced from the luminosity 
center. For concentrated clumps, the ratio is insensitive to the displacement. We note that 
we have not taken account into the extended nature of the galaxy. Had we compared the 
off-center clump's scattered intensity to the scattered intensity of an extended galaxy instead 
of a point source then the ratio would have been larger, c.f. §IIIa. 

e. Multiple Point Luminosity Sources 

In this section we examine the polarization pattern produced by the combination of two 
point sources embedded in a scale-free power-law scattering distribution (a = 1 in eq. (6)). 
For illustrative purposes, we take the two sources to be equal in luminosity and equidistant 
from the gas center as seen in projection. One source remains fixed in the plane of the sky 
that passes through the gas center, and the other source is placed at predetermined intervals 
perpendicular to this plane. In this example, the sources are oriented so that one lies at an 
angle of n/2 and the other at an angle of with respect to the gas center. 

The total polarization at any given point can be calculated numerically using eqs. (4), 
(9), and (10). First, the integrals in equations (4) are numerically evaluated for a grid 
of 100 x 100 lines of sight for a single source. Then, the contribution of a second source 
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equivalent to the first but rotated by n/2 is calculated and combined pointwise on the grid. 
Figure 4a shows the fractional polarization, the angle of polarization, and the intensity of 
scattered light for a configuration with both point sources lying in the plane of the gas 
center. As expected, the intensity has local maxima at the sources and at the center of the 
scattering distribution. The new feature, however, is that there are two regions of destructive 
interference seen as lobes of low polarization. It is easy to understand their cause. Consider 
the right triangle defined by the line connecting the two sources and the angles 7r/4 at 
each source. Clearly, near the triangle's third vertex the polarization angle of one source 
is perpendicular to the other and the flux from the sources is identical, so the scattered 
polarized intensity is zero. On the other hand, the largest polarization must occur midway 
on the straight line connecting the two sources and, also, at large distances from both. In 
both of these regions the polarization adds constructively and approaches the asymptotic 
value of 50% for a = 1. We conclude that multiple sources complicate the polarization 
pattern where the scattering intensities are roughly equivalent. 

The remaining plot, Figure 4b, shows the polarization pattern when the source at n/2 
(the "secondary source") is moved along the line of sight to a distance four times its projected 
distance from the gas center. As the source is moved out of the plane of the sky the scattering 
properties undergo noticeable changes. (It does not matter, however, whether the source is 
moved into the foreground or the background; both directions will show the same change.) 
First, the intensity of the scattered light around the secondary source decreases from its 
maximum value in the plane (which is consistent with the previous example of a single 
source). Second, the fractional polarization in the vicinity of the secondary source drops in 
this example from roughly 35% to 20% as line of sight distance goes to its maximum value 
(which is also consistent with our previous example of a single source). Finally, the area 
of interference of the two sources (the two low-polarization lobes) approaches the secondary 
source as it moves out of the plane. The lobes occur where both galaxies provide the same 
flux so that the parallel and perpendicular components of the light scattered by both sources 
negate one another. That projected location must be closer to the secondary source where 
it is moved out of the plane of the sky. 
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f. Summary 

Our examples show that the geometry of the source and of the electron distribution 
governs the polarization patterns in a simply understood fashion. One can envision many 
potential applications. It may be possible to constrain the three-dimensional position of 
luminosity sources by studying such patterns. For example, consider a luminous galaxy lo- 
cated off-center in a spherical-cluster potential that is filled with ionized gas. The galaxy's 
offset position along the line of sight with respect to the cluster center can be determined 
if enough of the light scattered by the cluster gas can be detected and mapped. Another 
potential application is to use scattering measurements to locate, identify and map inhomo- 
geneities, e.g. as gas becomes thermally unstable in a cooling flow. Finally, in a cluster with 
many galaxies present, the constructive and destructive effects in the polarized scattering 
lead to complicated neighbor-neighbor interactions. Detailed polarization maps may be used 
to model the three-dimensional distribution of galaxies in a cluster. A prime candidate is the 
cluster 2A 0335+096 (Sarazin, O'Connell, and McNamara 1992). This cluster has a central 
D galaxy with a bright companion nucleus at a radius of 6 kpc. Added to this are two 
other galaxies located at distances of 40 kpc. Throughout this region, ROSAT observations 
indicate that the gas density is slightly less than 0.1cm -3 . In the section that follows we 
consider additional physical effects that complicate the use of polarization as a potential 
observational tool. 

III. Estimates of Scattering From Actual Sources 

a. Extended Luminosity Sources 

We have previously considered the intensity of the polarized scattering, the angle of 
polarization, and the fractional polarization produced by point-like luminosity sources. In 
practice the luminosity source's spatial extent may also be important, e.g. if galactic emission 
is scattered by gas within a few core radii of the galaxy itself. In this section we consider the 
continuum limit for the luminosity source and characterize the differences in the scattered 
emission for a point source and an extended source. 
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Let pi be the luminosity density and n e be the scattering density, each assumed to be 
spherically symmetric about the origin. Without loss of generality, choose the line of sight 
in the z direction and py as the impact parameter. Let / be the distance along the line of 
sight and s the distance between (x,y,z) and (0,p, I). The scattered intensity and polarized 
intensity are then given by integrals over the volume and along the line of sight: 

I 
P 

where P = ^Q 2 + U 2 . 

We assume a spatial luminosity density of the form 

PL(r)= PL (0)(l+r 2 /r 2 c f , (28) 

where Pl(0) is the central luminosity density, and r c is the core radius. As a specific example, 
we choose (3 = —3/2, which implies that the surface brightness scales like r~ 2 power law 
for large radii. Such a power law is commonly observed for giant elliptical and cD galaxies 
(Oemler and Tinsley 1979). The total luminosity is dictated by the outer cutoff, which we 
take to be 100 r c . However, all the scattering results we discuss are essentially independent 
of the cutoff because they depend upon integrals of functions which fall much more steeply 
with r than Pi(r). We assume a radial distribution of the scattering medium given by 

n e (r) = n c (l + r 2 /r 2 c y . (29) 

We are primarily interested in the scattering in the halo; we include the core r c to prevent 
an unrealistic divergence at the origin (as occurs in eq. 6) As a specific example, we choose 
7 = —1/2, which is roughly appropriate when the scattering medium is part of a cooling 
flow (White and Sarazin 1987, Meiksin 1990), based on model fits to X-ray observations. 
There are complications in using these results to infer the scattering profile. In the models 
the X-ray intensity constrains (n 2 ) but (n e ) can be varied by adjusting the gas clumping. 
Since the scattered intensity is proportional to (n e ) this yields a range in J, which Wise and 
Sarazin (1992) have shown for the galaxy M87 can be as much as 40%. In addition, the loss 
of ionized gas by recombination and star formation alters the inferred scattering density. 
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Substituting eqs. (28) and (29) into eq. (27) we find 



I 
P 



3crj-p^(0)n c r^ r dxdydzdl 



647T 2 / 5 4( 1+5 2 + £2 + 5 2)/3( 1 + /2 + p2) 7 

where the dimensionless quantities x, y, z, I, s, and ft are x/r c , y/r c , etc., respectively. The 
integral has been evaluated using Romberg integration at 20 impact parameters. 

The calculated polarization as a function of r, II, is shown in Figure 5. Because of the 
four- dimensional nature of the integral the numerical computation is quite time consuming. 
Near the core radius up to 10 11 function evaluations per impact parameter were done. The 
relative accuracy was roughly a few percent in the inner regions. Note that II — > as r — > 0, 
since the source illuminates the line of sight symmetrically, and that II^l/2asr^oo 
(beyond the cutoff) when the point-source approximation is good. The figure shows that it 
is necessary to be many core radii (r > 10 r c ) from the center before II has risen to even 
half its asymptotic value. Figure 6 shows the ratio of / and P for the extended source to / 
and P for a point source. It can be seen that / and especially P are reduced significantly 
as the impact parameter approaches the core of the galaxy. Consequently, the intensity of 
polarized emission is much less than would be predicted given / and the asymptotic value 

of n. 

Since the quadratures are time consuming, we give two empirical rules for the total and 
polarized intensity as a function of radius for extended sources, 

J(r)^/ pt (r) L ^ (r)+ r L ^ (rc) (31) 



~s 2 + (l-~z) 2 

(y — P) 2 — % 2 



(30) 



and 



P(r) - lUr)^^ • (32) 



Here L cy i(r) and L sp /j(r) represent integrations of equation (28) giving the total luminosity 
within a cylinder and sphere of projected radius r, respectively. I p t(j') is the point source 
intensity (c./. eq. (7)) given by 

La T r c n c 

W = 87r 2 r 2 • ( 33 ) 
Our rules of thumb given above are accurate to within 30% for 1 < r/r c < 100. Inside the 
core of the galaxy the values of / and P change too rapidly to be represented accurately by 
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empirical laws. The total luminosity of the galaxy, idealized as a point source, provides a 
good fit only at large radii. For an extended source both the scattered intensity and polarized 
intensity are reduced by a significant fraction within about 10 core radii. 



In this section we predict the intensity of scattered light from cD galaxies. Such galaxies 
have two desirable properties: (i) they are embedded in massive cooling flows; and (ii) they 
are usually the most luminous galaxy in the cluster in which they reside. A luminous source 
in a dense scattering gas gives us the most favorable conditions for observing scattered light. 
We examine two galaxies in detail: M87 and NGC 1275. Each is associated with a cooling 
flow and each is the dominant galaxy in its cluster. We assume that the scattering gas is 
centered on the galaxy and that the galaxy acts as an extended source of luminosity. We 
make use of our results from §IIIa. 

NGC 1275 has an enormous cooling flow, ~ 300 M yr _1 , and the total mass of gas lying 
within 200 kpc is estimated to be 10 13 M (Fabian et al. 1981). We assume the luminosity 
density and electron scattering gas are given by eqs. (28) and (29), respectively. We take 
f3 = —3/2 and 7 = — 1/2 as discussed in §IIIa. Integrating eq. (29) and expressing n c in 
terms of the total mass, M to t within the radius r g , we find 



where m p is the proton mass. The 7/8 that appears in eq. (34) assumes the gas is composed 
of 75% ionized hydrogen and 25% doubly ionized helium by mass. The core radii for both 
scattering and luminosity densities are assumed to be 3.2 kpc, based on the V-band obser- 
vations of Oemler (1976). We find n c ~ 0.45 cm -3 . Ly is known from the absolute visual 
magnitude of -24.27 (Oemler 1976). Finally, we assume an outer cutoff in the luminosity 
density at r cut = 100r c and solve for p_c y (0). By integrating eq. (28) we find 



b. cD Galaxies 




(34) 




(35) 
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and determine PL y (0) = 2.96 x 10 23 ergs 1 cm 3 . Knowing PL y (0) and n c , it is now possible 
to calculate the total and polarized intensities for our model galaxy from eq. (30). 

Figure 7a shows the total and polarized intensities, along with 1% of the background 
intensity in V. The background is composed of galaxy plus sky emission and is a function 
of r. We include it as an effective lower bound that the total and polarized intensities must 
exceed for any hope of detection. As a practical matter, current instruments can measure 
an intensity > 1% of the background (Miller and Goodrich 1988). We take the sky intensity 
in V to be 21.9 mag arcsec" 2 (Pilachowski et al. 1989). The galaxy's contribution to the 
background decreases at shorter optical wavelengths, but the intensity of the scattered light 
is similarly diminished. As we will show in §IIIc, when the scattered light derives from a 
bluer source (such as an AGN) it is advantageous to observe in U or B. In this example, 
however, little can be gained by doing so. An additional difficulty in making measurements of 
polarization fractions below 1% of the galaxy emission is the grain-induced polarization that 
affects light traveling through the Galaxy's interstellar medium. Studying a large sample of 
nearby stars, Serkowski, Mathewson, and Ford (1975) found an empirical expression for the 
maximum polarization, II < 0.09 E(B — V), where E(B — V) is the differential extinction. 
For \b\ < 20° the range of E(B - V) is roughly 0.01 to 0.1 (Burstein and Heiles 1982), with 
the lowest values found toward the poles. Hence, the grain-induced polarization should be 
less than 1% of the galaxy's emission. Our results show that detecting the scattered light 
in NGC 1275 is currently impossible at optical wavelengths. At no point is / above the 
1% background level. Furthermore, the polarized intensity is a very small fraction of the 
total scattered intensity. The maximum of polarized emission occurs just beyond the core 
radius, only 27.8 mag arcsec -2 , or roughly 3 magnitudes fainter than the 1% level of the 
background. 

We include a second example with a different core radius and scattering gas density to 
suggest that it will generally be impossible to observe scattered light in the vicinity of a 
cD galaxy. This example, M87, has a maximum mass inflow rate of 4OM yr _1 . Here n c is 
based on the cooling flow model fit to Einstein X-ray data (White and Sarazin 1988). Oemler 
(1976) gives the core radius of 2.2 kpc and we match the number density of the cooling flow 
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model there. The result is n c = 0.75 cm~ 3 . The luminosity density is calculated in the same 
manner as above using eq. (35), for an absolute visual magnitude of -23.53 (Oemler 1976). 
We find pLy(O) = 4.6 x 10~ 23 ergs _1 cm~ 3 . 

Figure 7b shows the values of 1% of the galaxy plus sky intensity, J, and P for M87. 
Again, the scattered and polarized intensities lie below the 1% level. We conclude that the 
intensity of visible scattered light from a cD galaxy immersed in a cooling flow is too small 
to be observed. 

c. A cD Galaxy With an Isotropic Active Nucleus 

In this section we add an active nucleus to the cD galaxy. In our model we assume (1) 
that the luminosity of the central nucleus is equal to that of the surrounding galaxy in V; (2) 
that the emitted light from the nucleus is isotropic; (3) that the nucleus is a point source; 
(4) that the nuclear luminosity is constant; and (5) that emission is unpolarized. In the next 
section we consider an anisotropic flux from the central source. For relaxation of some of 
these restrictions we refer the reader to Wise and Sarazin (1992). 

In this section we will compare the scattering in U, B, and V bands. Thus, we assume 
(6) the parent galaxy has a spectrum of the form, 

vgal „ j is' 2 ' 1 , for v < 10 1483 Hz 
r v °C \ ,,-6.0 f_ „ ^ inl4.83i 



i/ _b - u , for v > 10 i4 83 Hz , 



and the active nucleus of the form 



I/" 1 , for v < 10 14 - 85 Hz 
z/°, for v > 10 14 - 85 Hz . 



These spectral distributions are rough fits to the "standard galaxy" taken from Malkan 
and Oke (1983) and to AGN's taken from Francis (1993). The U, B, and V intensities are 
determined by integrating the above distributions over the appropriate frequency ranges. 

We choose the same parameters as those for NGC 1275 except that we include a central 
point source of equal luminosity, L*y U = 1.54 x 10 45 erg s _1 . Far outside the core of the 
electron scattering distribution, eq. (7) implies that the scattered intensity from the active 
nucleus goes as 1/r 2 . Figures 8 show the total and polarized intensities for V, B, and 
U. For U and B we take the sky intensity to be 22.2 and 23.0 mag arcsec -2 , respectively 
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(Neizvestnyi 1982; Pilachowski et al. 1989). In V the total and polarized intensities exceed 
the 1% background level for radii less than 15 kpc and 4 kpc, respectively. Moving to U these 
radii become 19 kpc and 13 kpc, respectively. Within 10 kpc the measured polarization, P, 
exceeds several percent in U. This change is largely due to the decreased galaxy background. 

For all three bands, the detectable effect is due completely to the presence of the active 
nucleus. Obviously, if the nucleus were brighter than the value we assigned it, the values 
of / and P would be proportionally larger. Active galaxies can have luminosities up to 
and exceeding 10 47 erg s _1 (Weedman 1986), roughly 100 times the luminosity of the most 
luminous cD galaxy. At a luminosity of 10 47 erg s _1 the scattered / and P, for parameters 
similar to NGC 1275, could be detected out to roughly 200 kpc. Although a correlation 
exists between the luminosities of the active nucleus and the parent galaxy, we note that we 
have already adopted a near-maximum estimate for the luminosity of the cD (Hutchings, 
Crampton, and Campbell 1984) whose background serves to obscure the scattered light. 
Many more galaxies may be visible in scattered light; this depends upon how many have 
nuclei that outshine the rest of the galaxy and, at the same time, have sufficiently high 
Thomson optical depths. 



There have been recent suggestions that the brightest active galactic nuclei are those 
that are being beamed toward us (Barthel 1989), though there is some evidence against 
this suggestion (Boroson 1992). With this in mind we now investigate a model that has 
its radiation beamed at an angle relative to the observer's line of sight. Let L& be the 
apparent luminosity for an observer when the beam is aligned with the line of sight. Relative 
to the isotropic case, with L = L^, the intensity of scattered light will be diminished since 
the emissivity is confined to the cone of radiation. From eqs. (4) approximate expressions 
for the intensity of scattered and polarized light are 



d. Active Nucleus with Beamed Radiation 




^A0 sin 0(1 + cos 2 0) 



(36) 



and 




2 A6 sin 0(1 — cos 2 0) 



(37) 
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Here p is the impact parameter, A8 is the angular width of the beam, and Iq = 3Li,n c aTr c 
/6Att 2 . These formulae are accurate when the beam is not pointed near the line of sight and 
A9 is small. As an example assume the gas distribution of NGC 1275, = 10 47 erg s _1 and 
A9 = 15°. If the beam is perpendicular to the line of sight, / ~ 9.0 x 10~ 13 /p 2 (kpc) erg s _1 
cm~ 2 arcsec~ 2 or 18.4 + 51ogio[p(kpc)] mag arcsec -2 and II = 1.0. Thus the scattered ra- 
diation from the beam could be observed out to ~80 kpc before it drops below the 1% 
background level. Even if the beam is not pointed toward us, its scattered light may be 
visible and may be used to determine the orientation of the beam (Fabian 1989; Wise and 
Sarazin 1992). 

e. Spectral Lines 

The presence of strong spectral lines from an AGN provides another possible strategy for 
detecting scattered light. If the brightness of the emission line from the AGN is significantly 
larger than the brightness of the continuum of the AGN plus galaxy background, then it 
may be advantageous to use a filter narrower than U, B, or V centered on the line frequency 
to enhance the signal-to-noise ratio of an observation. Scattered light from regions far from 
the nucleus could be more easily seen because background contamination would be reduced. 

There are several complications. First, the scattered line spreads in frequency because 
the photons scatter off a hot gas, c.f. § Ha. The expected temperatures in cooling flows 
range from 10 6 K (inner few kpc) to 10 8 (outermost regions) (Fabian et al. 1991) so that a 
scattered optical line with a wavelength of ~ 5000A will be smeared from 60A to 600A. If 
the nuclear continuum is flat over similar bandwidths, the smearing reduces the brightness 
ratio of the nuclear emission line to the nuclear continuum. 

Observations of AGN's show that the strongest lines are roughly 3 times the AGN 
continuum. Many of these lines such as Ha have velocity widths of a few thousand km s^ 1 or 
~ few xlOOA (Osterbrock and Mathews 1986). Narrower lines such as [OIII] are less useful: 
their maximum intensity per frequency interval is comparable to Ha, but the smearing will 
reduce their intensities relative to the continuum by a much larger factor. 

To make an estimate we proceed as follows. Some emission lines such as Ha arise almost 
solely from the photoionization of clouds in the inner 100 pc of the nucleus. In section 
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§IIIc we have shown that the most promising region for detecting scattered light is within 
the inner 20 kpc of the nucleus. At these radii the gas temperature will be ~ 10 6 K and 
Ha will be smeared ~ lOOA. Since the smearing is less than the intrinsic line width, the 
ratio of scattered line to scattered nuclear continuum will be roughly 3. If the background 
from the galaxy is roughly flat over the same bandwidth, then the filtered intensity may be 
enhanced by up to a factor of 3 compared to the broad-band case. Thus, strong emission 
lines emanating from an AGN could be marginally helpful for observing scattered light in 
regions of a cooling flow outside of the nucleus. 

One complication is that many cooling flows exhibit the same optical emission lines 
as found in AGN's (Stockton and MacKenty 1987). In this case, the galactic background 
includes a component that is not reduced by using a filter and the usefulness of the emission 
lines is decreased. Still, one could distinguish the origin of this emission by polarization 
measurements and possibly by subtracting the narrower background emission line of the 
cooling flow. 

f. OS-Center Gas Distribution 

We have shown that scattered and/or polarized light can currently be detected in a cD 
galaxy with an active nucleus of luminosity roughly equal to or greater than that of the rest 
of the galaxy. In this section we use the results from §IId to examine the situation when part 
of the scattering gas is not centered on the cD galaxy. For a galaxy of fixed luminosity, a 
scattering halo typical of a cD (a = 1) and a clump of fixed column density we showed that 
the ratio of the clump scattered brightness to that of the halo increases with distance from 
the galaxy center. Two additional factors must be taken into account to assess the visibility 
of scattering from such off-center clumps: (i) the diffuse background of the galaxy and (ii) 
the intrinsic sky level. 

A possible example of an off-center density distribution is the radio source PKS 0745- 
191 (Fabian et al. 1985). This galaxy is immersed in a hot gas that is cooling at a rate of 
up to 1000 M Q yr _1 . The X-ray image of this source extends to roughly 200 kpc with a 
peak of emission that lies 30 kpc from the center of the galaxy (indicating a possible density 
maxima). From the optical observations of Fabian et al. (1985) we arrive at the following 
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values for the galaxy: r 9 c al = 9 kpc, r cut = 45 kpc, my = 16, and Ly = 1.09 x 10 45 erg s _1 . 
The cutoff radius of the galaxy is uncertain but its actual value does not change our results 
significantly. From X-ray data the average density within 30 kpc of the intensity peak is 
0.06 cm~ 3 (Fabian et al. 1985). We take the scaling radius of the off-center density peak 
to be equivalent to the core radius of the cD galaxy, rf = 9 kpc. Then, assuming for the 
scattering gas a power law identical to eq. (6) where a — 1, we arrive at = 0.13 cm -3 . 
The luminosity is assumed to be emanating from a AGN for the purpose of calculating the 
scattering. (Consulting Figure 6, we expect the errors in the total scattered intensity to 
be less than a factor of two.) For the purpose of determining the background emission we 
assume the galaxy's luminosity density is extended, following eq. (28) where f3 — — 3/2 . 

Using the information from the previous paragraph we now estimate the intensity a 
distance rf from the scattering gas center. We will assume that the gas center is coplanar 
with the galaxy center and the clump gas has a — 1. We find the intensity of the polarized 
emission, P, to be 27.5 magnitudes per square arcsecond, about one magnitude brighter 
than the examples in §IIIb at the equivalent radius. In our power-law model the intensity 
diverges as the density singularity is approached (c.f. §IId); the actual maximum is dictated 
by the clump's core radius which, of course, is unknown. Overall, there is a region roughly 
10 kpc in diameter near the density peak with I and P above the 1% background level. 
Thus, we can only conclude that if the core radius < 10 kpc scattering may be visible. If 
the gas center lies out of the plane of the galaxy then the polarized intensity will be reduced 
significantly relative to a gas center in the plane. Thus, it may be possible to detect relatively 
dense regions of hot gas residing near a luminous galaxy. If detected, it may be possible to 
determine the position along the line of sight from the fractional polarization. 

Finally, we point out that the presence of an active nucleus with a luminosity greater than 
that assumed above lifts the level of scattered light proportionally. Scattering may be used 
to probe the galaxy's gas and cooling flow and so represents a potentially complementary 
source of information to that provided by X-ray observations. In a cooling flow, for example, 
one would like to know where the infalling gas is converted into stars plus the length scale 
and growth rate of thermal instability. These angular scales are typically smaller than the 
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resolution of X-ray telescopes, but easily within the range of optical instruments. 

IV. Conclusions 

In this paper we have outlined the nature of scattering in a hot gas surrounding a 
luminous galaxy or an active galactic nucleus for a variety of simple models. We began 
by considering a point-like luminosity source embedded in scale-free electron power laws. 
We derived asymptotic estimates for the scattering intensity and polarization fraction in 
the neighborhood of the density peak and the luminosity source. We complemented these 
analytic results with numerical calculations of the scattering intensity, polarization fraction, 
and angle of polarization for single and multiple point-like luminosity sources. The large-scale 
patterns illustrate the effects of interference between sources. For an extended luminosity 
source we showed that the interference decreased the polarization fraction in the vicinity of 
the source. 

We applied our results to investigate several potential types of scattering in the vicinity 
of galaxies. For a typical cD (without an active nuclei) the scattered galactic emission is 
too weak to be detected with current instruments at optical wavelengths. However, those 
galaxies with an active nucleus and/or off-center gas are possibilities. 

An active nucleus increases the flux of scattered radiation and is so small that interference 
effects are negligible. We show that when a galaxy has an active nucleus equal in luminosity 
to the parent galaxy, the total and polarized scattered intensities in the inner 5-30 kpc 
exceed the 1% background level. This is a conservative estimate since the luminosity of 
many nuclei exceed their parent galaxies by factors of 10 to 100; for these more luminous 
examples the scattered and polarized light becomes visible on much larger size scales. For 
typical galaxy and AGN spectra, it is significantly easier to measure the scattering at short 
wavelengths (i.e. U). If an active galaxy has an obscuring dust torus or is beaming its 
radiation in a direction away from the line of sight, it may be possible to determine the 
orientation of the torus/beam by searching for the cone of scattered light. There is only 
modest advantage to using a narrow-band filter centered on an emission line because the 
typical continuum associated with the active nucleus is only a factor of 3 less in intensity 
than the scatter-broadened emission lines. 
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For galaxies with off-center gas distributions, detection is easier because of the reduced 
background from the galaxy at the gas center and a high degree of polarization. The polar- 
ization fraction reaches a local maximum as the impact parameter approaches the off-center 
distribution. In our example, PKS 0745-191, the scattered and polarized intensity were 
greater than 1% of the background intensity of the galaxy and sky for a region roughly 
10 kpc in diameter. When an active nucleus is present both the intensity of the scattered 
radiation and the size of the region grow. Hence, clumps of ionized gas in the environment 
surrounding the galaxy should be detectable by searching for regions of polarized emission. 

We conclude that there are a number of promising observational possibilities offered by 
polarized scattering that occurs near galaxies, especially those with AGN's. 
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Figure Captions 

Figure 1. A photon leaves the source in the direction k, is scattered through an angle 9 into 
the line-of-sight h. The polarization bases of the the unscattered light, e\ and e2 (see 
§IIa), and the scattered light, ej_ and ey, are illustrated. 

Figure 2. The parallel and perpendicular bases for the scattered light are displayed as seen 
by the observer. 

Figure 3a. Polarization fraction (length of the line segments), the angle of polarization 
(orientation of the line segments), and the intensity of the scattered light (contours) for 
a luminosity source coplanar with the gas center. Here the luminosity source is coplanar 
with and lr c to the right of the gas-density origin. The contour interval is one-half 
magnitude (10 2 ). 

Figure 3b. Contours of polarization fraction are shown for the same geometry as in Figure 
3a. The contour interval is .04. 

Figure 3c. Same as 3a but the luminosity source is 4r c out of the plane of the gas-density 
origin. 

Figure 4a. Same as 3a but with two luminosity sources coplanar with the gas-density origin. 

Figure 4b. Same as 3a but with two luminosity sources, one coplanar with the gas-density 
origin and one 4r c out of the plane. 

Figure 5. Polarization as a function of radius for an extended source given by eq (30). 

Figure 6. The ratios of the extended to point source scattered intensities for / and P. 

Figure 7a. The scattered brightness profile of NGC 1275 (I and P) is compared to the 
threshold for detection, which is 1% of the sum of the extended diffuse galactic emis- 
sion and the sky. /, P, and the 1% background level are represented by the dashed, 
dotted, and solid lines, respectively. At no point is the scattering detectable in optical 
wavelengths with current instruments (see §IIIb) . 

Figure 7b. Same as 7a for M87. 

Figure 8. Same as 7a except that a central point source of luminosity Ly = 1.54 x 10 45 ergs _1 
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(equal to that of the background galaxy) is included in the calculation for three wave- 
length bands, U, B, and V. /, P, and the 1% background level are represented by the 
dashed, dotted, and solid lines, respectively. Note that / and P exceed the detection 
threshold for this case. 
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